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Abstract: Net ecosystem exchange (NEE) of carbon dioxide (CO2) was measured in a cool 
temperate peatland in northwestern Turkey on a continuous basis using eddy covariance 
(EC) sensors and multiple (non-)linear regression-M(N)LR-models. Our results showed 
that hourly NEE varied between -1.26 and 1.06 mg CO2 m~ 2 s 1 , with a mean value 
of 0.11 mg CO2 m~ 2 s~ . Nighttime ecosystem respiration (R E ) was on average measured 
as 0.23 ± 0.09 mg C0 2 m" 2 s" 1 . Two best-fit M(N)LR models estimated daytime R E 
as 0.64 ± 0.31 and 0.24 ± 0.05 mg C0 2 m" 2 s" 1 . Total R E as the sum of nighttime and 
daytime R E ranged from 0.47 to 0.87 mg CO2 m s , thus yielding estimates of gross 
primary productivity (GPP) at -0.35 ± 0.18 and -0.74 ± 0.43 mg C0 2 m" 2 s~\ Use of EC 
sensors and M(N)LR models is one of the most direct ways to quantify turbulent CO2 
exchanges among the soil, vegetation and atmosphere within the atmospheric boundary 
layer, as well as source and sink behaviors of ecosystems. 
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1. Introduction 

Though spatially small (5% of the terrestrial biosphere) compared with most other ecosystems [1,2], 
peatlands play a significant role in carbon (C) and water metabolism of the World. Understanding and 
quantifying C dynamics of peatlands are crucial to prediction of responses to global climate change and 
rehabilitation of peatlands under the increasing magnitude and rate of human-induced disturbances. Eddy 
Covariance (EC) sensors are one of the most direct ways to measure and estimate turbulent carbon 
dioxide (CO2), water vapor and energy fluxes exchanged among soil, vegetation, and atmosphere within 
the atmospheric boundary layer. The use of the EC method and sensors on a long-term and continuous 
basis across the World has led to the establishment of an integrated global network for standardization of 
flux tower activities (called FLUXNET) and a network for standardization and development of spectral 
sensors toward bridging the gap between remote and proximal sensing (called SpecNET). 

Our study area in Turkey, the Yenicaga peatland area that at one time occupied 240 km , has 
diminished to less than 30 km due to drainage, cultivation, afforestation, and peat mining [3,4]. There 
is a lack of information about C metabolism of peatland ecosystems in Turkey, and this study is the 
first to comprehensively determine C dynamics and components in one of the remaining major 
peatlands. The objective of this study was to quantify the rate, magnitude, and timing of CO2 exchange 
between the atmosphere and Yenicaga peatland using the EC sensors and multiple (non-)linear 
regression-M(N)LR-models. 

2. Study Site and Methodology 

2.1. Description of Study Site 

The Yenicaga peatland is located about 38 km east of the city of Bolu (40°47N', 32°1E) in the 
western Black Sea region of Turkey (Figure 1). 

Figure 1. Location of the study site "Yenicaga peatland" in northwestern Turkey. 




2 

An EC flux tower site was installed about 1 km north of Lake Yenicaga (18 km ) at the elevation 
of 988 m above seal level on July 12, 2010. The climate in the Yenicaga region is classified as cool 
temperate, with a mean annual temperature and precipitation of 10.2 °C and 538 mm, respectively [5]. 
The Yenicaga peatland is reported to contain Devonian and cretaceous limestone, basaltic tuff, lava, 
and olistolites, with the uppermost layer consisting of tertiary and quaternary formations [5]. The 
natural vegetation types of the Yenicaga peatland consist of the following dominant communities: 
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(1) Phragmites australis and Typha domingensis (41.2 ha); and (2) Ranunculus lingua, Acorus 
calamus, Najas marina, Pedicularis palustris, Senecio paludosus (107.8 ha) [6]. The mean vegetation 
height around the flux tower is about 0.5 m, and the terrain observed around the flux tower exhibits 
flat and uniform grasslands. 

2.2. Eddy Covariance Flux and Ancillary Measurements 

Net ecosystem exchange (NEE) rates of carbon dioxide (F c , mg m s ) in Yenicaga peatland were 
estimated using an eddy covariance (EC) system consisting of an open-path CO2/H2O gas analyzer 
(LI-7500, Licor Inc., Lincoln, NB, USA), a 3-D sonic anemometer/thermometer (CSAT3, Campbell 
Scientific Inc., Logan, UT, USA), a data logger (CR3000, Campbell Scientific Inc.), and 
a 3-m tower on which EC flux sensors were mounted. The distance between the LI-7500 and CSAT3 
sensors was 0.15 m, with CSAT3 oriented towards the prevailing wind direction (an azimuth angle 
of 30° from true north) and LI-7500 vertically rotated 15° towards the footprint. 

Eddy fluxes and associated signals were recorded at 10 Hz, block averaged over one hour (h) and 
corrected for the effects of fluctuations in air density on CO2/H2O fluxes (F c wp i and LE wp i with WPL 
correction) through the online flux computation. EC data were collected swapping two 2-GB Compact 
Flash cards at 14-to- 18-day intervals. The net radiation (R n ), downwelling and upwelling longwave (4 
to 50 jum) and shortwave radiation (0.2 to 4 /mi) (i?i_dn, Ri_ up , R s dn, and R s _ up , W m~ 2 ) were measured 
using Kipp & Zonen CNR-4 net radiometers (Kipp & Zonen USA Inc., Bohemia, NY, USA). Air 
temperature (r a , °C), and relative humidity (RH, %) were sampled using HMP45C probe (Vaisala, 
Finland). Precipitation (PPT, mm), evapotranspiration (ET, mm), soil water content (SWC, %), and 
mean, maximum and minimum soil temperature (ST, ST max , and ST m i n , °C) were measured on a hourly 
basis using ET107 weather monitoring station (Campbell Scientific Inc.). Photosynthetically active 
radiation (PAR) was estimated from net shortwave radiation (R S D ) using a conversion factor of PAR: 
R s n = 0.5 [7-9]. The values of CO2 fluxes in unit of mg m s were also converted to unit of 
kg C ha _1 day _1 , based on the following conversion ratios of mg:kg = 10 6 ; s:day = 86,400; 
m 2 :ha = 10,000; and C0 2 :C = 44/12. 

2.3. Data Processing and Analyses 

As with the conventional meteorological sign notation, downward and upward fluxes are 
considered negative (— ) and positive (±), respectively. The first step of data processing involved the 
removal of erroneous spikes and their associated CO2 fluxes when latent heat flux (LEwpi) < -100 
or > 800 W m~ 2 ; sensible heat flux (H s ) < -150 or > 500 W m~ 2 ; and precipitation events occurred [10]. 
For both hourly daytime and nighttime CO2 fluxes, descriptive statistics were given, and best-fit 
cumulative distribution function (CDF) was selected. Tukey's multiple comparison was performed 
after one-way analysis of variance (ANOVA) for the entire dataset in order to test significant 
differences in hourly, nighttime versus daytime and monthly means. All the statistical analyses were 
performed with Minitab 15.1 (Minitab Inc. 2006). 

Second, EC data were separated into daytime (R a and/or R & dn > 10 W m ) and nighttime (R n and/or 
R SJ in < 10 W m ) periods since EC data are more reliable during daytime hours than during nighttime 
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hours, and nighttime EC data ( F" 1 ^ ) can be used to estimate daytime as well as nighttime ecosystem 
respiration (both plant and soil respiration) (R E = F" l ^ x ). Negative night CO2 flux data were deleted as 
no gross primary productivity (GPP = 0) occurs during the nighttime. The site-specific threshold value 
of friction velocity (u*) was determined as 0.03 m s -1 below which low vertical wind velocity led to 
underestimation of the nighttime CO2 fluxes [11]. Likewise, nighttime periods where horizontal wind 
velocity was less than 1 m s -1 were removed from the dataset. Finally, night fluxes where CO2 density 
data had the standard deviation of > 14 mg m (0.6 jomole m ) were eliminated from further 
analyses [12]. Multiple (non-)linear regression models were fitted to the resultant night CO2 dataset, 
and best-fit M(N)LR models with and without the inclusion of temporal variables (hour, and month) 
were chosen using best subsets procedure (low Mallows' C p , high adjusted R , and low SE). The 
following soil respiration equations of RothC [13] and CENTURY [14] models were also used to model 
daytime Re. 



/(*soii) = r 47 " 9 1()6 x for RothC model 

1 + exp 



(1) 



vST + 18.3 j 

f(R soil ) = 0.56 ± 0.465*arctan(0.097*(ST-15.7)) for CENTURY model ( 2 ) 

where 7? S0 ii is soil respiration (mg CO2 m s ), and ST is soil temperature (°C). Net ecosystem 
exchange (NEE, mg CO2 m s ) can be expressed as follows: 

NEE = -GPP ± R E (3) 

Hourly GPP values were estimated as the sum of NEE and daytime Re. Daytime Re was obtained 
extrapolating nighttime Re fluxes to the remaining daytime fluxes, based on M(N)LR models. 

3. Results and Discussion 

3.1. Diurnal CO 2 Fluxes 

Descriptive statistics of despiked and averaged EC data for the period of July 12 to October 17, 2010 
indicated that mean hourly CO2 fluxes above the canopy had a larger temporal variability (CV = 498) 
than the rest of the measured variables and varied between -1.5 and 1.5 mg m~ s~ (Table 1). During 
the study period, the study site cumulatively received 92 mm PPT and lost 305 mm water vapor via 
ET. The mean Bowen ratio of 7/ s :LE wp i quantifying the evaporative demand of an environment was 
estimated at 0.37, thus characterizing the Yenicaga peatland as a mesic environment. A logistic CDF 
with a location of 0.05197 and a scale of 0.1596 appeared to fit our daytime and nighttime F c wp i data 
fairly well as follows (r = 0.99; n = 1884; P < 0.001) (Figure 2): 

CDF logistic = l + e -(F c _ wpl -v)/s (4) 

where ju and s are location and scale parameters of the logistic CDF, respectively. The fitted logistic 
CDF can be used to estimate percentiles for the CO2 flux data (Figure 2). 
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Table 1. Descriptive statistics of eddy covariance and meteorological data collected on an 
hourly basis in the Yenicaga peatland. 



Variable 


n 


mean 




LV 


min 


median 


max 


^c_w P i (mg C0 2 m s ) 


1 O O A 

1,884 


0.06 


n 1 n 

0.30 


a no 

498 


1 en 

—1.50 


n nzi 

0.06 


1 c 1 

1.53 


H S (W m ) 


1 O O A 

1,884 


"> O Af\ 

38.40 


58.44 


152 


ho 0 £. 

-78.86 


i n 1 

3.01 


no on 

238.80 


T TJ SWT 

Lh wp i (W m ) 


1 OO A 

1,884 


102.69 


1 1 n 7o 

130. /y 


1 T7 

12 / 


ni on 
—92.89 


a 1 /in 
41.40 


51 /.29 


(m s ; 


1 OO A 

1,884 


n 1 0 
0. 18 


n 1 1 
(J. 12 


O8 


n n 1 
0.01 


n 1 /i 
0. 14 


n eh 
0.5 / 


CO2 density (mg m 3 ) 


1 OO A 

1,884 


633.34 


1 18.19 


19 


O /1 1 A E 

243.45 


cno EH 
598.5 / 


1 EHH H 

1,5 / 1 .15 


H 2 0 density (g m 3 ) 


1 OO A 

1,884 


12.34 


/I HH 

4. / / 


1 n 
39 


—2.62 


1 1 C\A 

12.94 


AH Ol 

4 /.82 


?a ( C) 


1 OO A 

1,884 


18.15 


t nc 

7.05 


i n 

39 


-0.67 


17.95 


1 /i nn 

34.00 


D TJ /O/ \ 

KH (%) 


1 OO A 

1,884 


H\ E A 

/1. 54 


22.30 


i 1 
3 1 


1 fin 
16.09 


H 1 1 

/6.13 


1 m on 
10 /.80 


WD (degree) 


1 OO A 

1,884 


219.91 


127.29 


c 0 

58 


n n a 

0.04 


^ £c ni 

265.93 


1 cn on 

359.89 


Wa(ms ) 


1 OO A 

1,884 


1 T5 

1. /3 


1 /I A 

1.40 


«1 


n n/i 
0.04 


1 1 n 
1.19 


6. 14 


«s„up(W m ) 


1 OO A 

1,884 


1 ^ on 

2o3.oU 


1 fin n 1 


11/ 


a n 1 
—4.91 


1 (\A 1 O 

104. 18 


no 1 0 0 
981.88 


« s _dn(Wm ) 


1 OO A 

1,054 


A A 1/1 

44. 14 


AO OO 

48.88 


111 
111 


n 1 n 

-0. 19 


18.63 


1 C 1 Cf 

153.55 


^(WnT 2 ) 


1,884 


386.85 


34.95 


9 


294.45 


386.94 


474.95 


^(Wnf 2 ) 


1,884 


413.84 


47.55 


11 


307.63 


405.96 


531.50 


^ n (W nf 2 ) 


1,884 


192.66 


245.44 


127 


-35.80 


68.69 


766.03 


ET (mm h" 1 ) 


1,680 


0.14 


0.21 


124 


0.0 


0.02 


0.75 


PPT (mm h" 1 ) 


1,680 


0.04 


0.36 


825 


0.0 


0.0 


7.8 


ST m i n ( C) 


1,680 


19.66 


3.92 


20 


8.11 


20.03 


26.18 


ST max ( C) 


1,680 


19.84 


3.96 


20 


8.19 


20.24 


27.24 


SWC (%) 


1,680 


75.37 


8.73 


12 


62.98 


73.75 


89.90 



F c W pi = WPL-corrected CO2 flux; H s = sensible heat flux; LEwpi = WPL-corrected latent heat flux; 
u* = friction velocity; T & = air temperature; RH = relative humidity; WD = wind direction; WS = wind 
speed; R s up = upwelling shortwave radiation; R s dn = downwelling shortwave radiation; 
^?i_up = upwelling longwave radiation; Ri dn = downwelling longwave radiation; R n = net radiation; 
ET = evapotranspiration; ST min = minimum soil temperature; ST max = maximum soil temperature; 
SWC = soil water content; CV = coefficient of variation; and SD = standard deviation. 

Figure 2. Logistic cumulative distribution function (CDF) fitted to both daytime and 
nighttime CO2 fluxes for the Yenicaga peatland: (Location (p) = 0.05197; 
scale (s) = 0.1596; r = 0.99; n = 1884; P < 0.001). 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



F c . wpl (mg C0 2 m 2 s ') 
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Quartic functions (an equation of fourth degree) fitted to both entire and monthly F c wp i datasets as a 
function of the explanatory temporal variable (h) provided a meaningful representation of CO2 fluxes, 

2 

with a R a dj value range of 70.9% in July to 3 1 .3% in October (Figures 3 to 7). 

Figure 3. A quartic function fitted to the entire F c _ wp i dataset in the Yenicaga peatland for 
the period of July 12 to October 17, 2010: F c wpl (mg C0 2 m" 2 s" 1 ) = 0.2807 h - 0.0602 h 2 
± 0.0038 h 3 - 7 x 10" 5 h 4 ; (i? 2 adj = 0.5368; SE: 0.2016; n = 1884; P < 0.001). 



2 T 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 



Local hour (h) 

Figure 4. A quartic function fitted to F c _ wp i dataset for Yenicaga peatland in July of 2010: 
F c wp i (mg C0 2 m" 2 s" 1 ) = 0.3575 h - 0.077 h 2 ± 0.0049 h 3 - 9 x 10" 5 h 4 ; {R 2 aA} = 0.7087; 
SE = 0.1905; n = 412; P < 0.001). 

1.5 T 




-1 -L 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 

Local hour (h) 
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Figure 5. A quartic function fitted to F c _ wp i dataset for Yenicaga peatland in August 
of 2010: F c wpl (mg C0 2 m~ 2 s" 1 ) = 0.3242 h - 0.0696 h 2 ± 0.0044 h 3 - 9 x 10" 5 h 4 ; 
(i? 2 adj = 0.5963; SE = 0.1957; n = 620; P < 0.001). 




Figure 6. A quartic function fitted to F c _ wp i dataset for Yenicaga peatland in September 
of 2010: F c _ wpI (mg C0 2 m" 2 s" 1 ) = 0.2297 h - 0.0488 h 2 ± 0.0031 h 3 - 6 x 10" 5 h 4 ; 
(i? 2 adj = 0.4726; SE = 0.1875; n = 603; P < 0.001). 



1.5 
1 

T ^ 0.5 

^ u 

£ -0.5 



-1.5 
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* ♦ V ft ft * ♦ 
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Figure 7. A quartic function fitted to F c wp i dataset for Yenicaga peatland in October 
of 2010: F c wpl (mg C0 2 m" 2 s" 1 ) = 0.1951 h 



0. 0423 h 2 ± 0.0027 h 3 



5 x 10" 5 h 4 ; 



adj 



0.3127; SE = 0.2074; n = 249; P < 0.001). 




+-4-t-+-*~*-! 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 

Local hour (h) 

3.2. Multitemporal Comparisons of CO 2 Fluxes 

The entire dataset of daytime and nighttime CO2 fluxes (n = 1,884) was used for a multiple 
comparison according to hours, nighttime versus daytime, and months, based on Tukey's test 
following one-way ANOVA (Table 2). A comparison of hourly mean F c _ wp i values showed that the 
Yenicaga peatland acted as a CO2 sink for the daytime periods between 9:00 AM and 17:00 PM and 
acted as a CO2 source for the periods between 18:00 PM and 8:00 AM. On average, the Yenicaga 
peatland had minimum and maximum CO2 fluxes ranging from -0.05 ± 0.21 mg CO2 m s at 9:00 AM 
to -0.28 ± 0.18 mg C0 2 m" 2 s" 1 at 12:00 PM as a C0 2 sink and from 0.05 ± 0.12 mg C0 2 m" 2 s" 1 
at 18:00 PM to 0.38 ± 0.29 mg C0 2 m" 2 s" 1 at 5:00 AM as a C0 2 source, respectively (Table 2). 

Mean daytime CO2 flux (-0.11 ± 0.22 mg CO2 m s ) significantly differed from mean nighttime 
C0 2 flux (0.26 ± 0.23 mg C0 2 m" 2 s" 1 ) (P < 0.001). All the monthly mean C0 2 fluxes were positive 
(upward into the atmosphere) and varied between 0.03 ± 0.35 mg CO2 m s in July 
and 0.07 ± 0.25 mg CO2 m s in September. However, the monthly mean CO2 fluxes did not appear 
to significantly differ from one another (Table 2). The best-fit MNLR model was derived from the 
entire dataset of nighttime and daytime F c wp i by exploring possible interaction effects among the EC 
and meteorological variables as follows: 



F c wp i (mg C0 2 m 2 s *) 



0.029586 h 2 ± 0.0018587 h 3 



0.001394 month ± 0.13776 h 
0.0000358 h 4 - 0.05387 daytime - 0.003694 H s ± 0.00001 183 H s 2 ± 0.016833 log(// s 2 ) 
0.0008682 LE wpl - 0.00000059 LE 2 ± 0.038155 log(LE 2 ) - 0.000013 P a - 0.009471 T a - 
0.000879 RH ± 0.00022167 RH*r a - 0.007483 R s dn ± 0.0000339 R s _ dn *Ri _u P - 
0.00000001 RH*T a *i? s dn *Ki up ± 0.009735 U z *R sdn 
(i? 2 adj = 62.6%; SE = 0.181; n = 1882; P < 0.001) 



(5) 
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where the variable "daytime" was used as an indicator variable coded as 1 and 0 for daytime and 
nighttime, respectively, and the sign "*" between the variable notations denotes two-to-four- way 
interactions. 



Table 2. A multiple comparison of mean F c _ wp i fluxes (mg CO2 m s ) in the Yenicaga 
peatland by different temporal scales (P < 0.001). 



Local hour 


n 


mean 


SD 


Comparison based on 95% confidence 
intervals 


1:00 


57 


0.254 


0.248 


("-*-) 
(--*") 
R") 
(-*-) 
(-*") 

(-*--) 

(-*--) 


2:00 


53 


0.313 


0.242 


3:00 


55 


0.337 


0.277 


4:00 
5:00 
6:00 


69 
65 
72 


0.320 
0.382 
0.276 


0.229 
0.296 
0.215 


7:00 


85 


0.209 


0.256 


8:00 


76 


0.087 


0.295 


(-*-) 


9:00 


82 


-0.053 


0.210 


(-*-) 


10:00 


90 


-0.171 


0.187 


(-*-) 

(-*-) 
(-*-) 

IJ 
(-*-) 
(-*-) 
(-*-) 

(-*-) 

(-*-) 


11:00 
12:00 


89 
91 


-0.242 
-0.280 


0.203 
0.187 


13:00 
14:00 
15:00 


91 
91 
89 


-0.237 
-0.191 
-0.200 


0.122 
0.168 
0.156 


16:00 


88 


-0.144 


0.158 


17:00 


90 


-0.069 


0.101 


18:00 


91 


0.052 


0.128 


(-*-) 


19:00 


89 


0.171 


0.145 


(-*-) 


20:00 
21:00 
22:00 


91 
85 
69 


0.222 
0.262 
0.275 


0.117 
0.149 
0.193 


(--*-) 
(-H 
(-*--) 


23:00 
24:00 


72 
54 


0.248 
0.267 


0.212 
0.266 


+ + + +-- 

-0.20 0.00 0.20 0.40 



Day vs. night 

Daytime 
Nighttime 

Month 

July 
August 
September 
October 



n 
1030 
852 

n 

412 
620 
603 

249 



mean 
-0.111 
0.266 

mean 
0.039 
0.054 
0.077 
0.064 



SD 
0.222 
0.235 

SD 
0.351 
0.307 
0.257 

0.248 



(*) 



(*) 



+— 
-0.12 



— +- 
0.00 
-* 



-+-- 
0.12 

-) 



-) 



0.24 



-) 



0.025 0.050 0.075 0.100 



3.3. Ecosystem Components of Carbon Metabolism 



Nighttime Re was modeled using (1) best-fit M(N)LR models with/without the forced inclusion of 
temporal variables (hour and month), and (2) ST-dependent 7? so ii equations [Equations (1) and (2)] of 
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RothC and CENTURY models with/without the forced addition of SWC. All the R E models were built 
with the intercept forced through zero as follows: 



(6) 



Re (mg C0 2 m" 2 s" 1 ) = 0.0004922 h - 0.021 1 1 month ± 0.13747 r sonic - 0.1496 T a ± 0.01 119 
ST min ± 0.0032329 SWC ± 0.5456 u* ± 0.5201 U z - 0.04903 WS - 0.00061 P a ± 0.0021203 R n 

-0.0027481 H s 
(R 2 ad] = 73.8%; SE = 0.0473; n = 203; P < 0.001) 

Re (mg C0 2 m" 2 s" 1 ) = 0.005383 r somc ± 0.0026921 SWC ± 1.0846 u* ± 0.0733 U z - 0.09595 

WS ± 0.0009581 H s - 0.0003665 RH 
(i? 2 adj = 70.8%; SE = 0.0677; n = 207; P < 0.001) 

Re (mg C0 2 m" 2 s" 1 ) = 0.33084 CENTURY 
(i? 2 adj = 42.0%; SE = 0.0708; n = 207; P < 0.001) 

Re (mg C0 2 m" 2 s" 1 ) = 0.085842 RothC 
(tf 2 adj = 44.3%; SE = 0.0686; n = 207; P < 0.001) 

Re (mg C0 2 m" 2 s" 1 ) = 0.071387 RothC ± 0.0005546 SWC 
(R 2 adi = 53.0%; SE = 0.0684; n = 207; P < 0.001) 

Re (mg C0 2 m" 2 s" 1 ) = 0.28185 CENTURY ± 0.0004839 SWC 
(i? 2 adj = 53.3%; SE = 0.0708; n = 207; P < 0.001) 

where T^mc is sonic temperature (°C), U z is vertical wind speed (m s" 1 ), and P a is air pressure (kPa). 
Nighttime Re Equations (6) and (11) were used to estimate daytime RE values which in turn led to the 
estimation of -GPP component from Equation (3) as can be seen in Figures 8 and 9. 



(V) 

(8) 
(9) 
(10) 

(11) 



Figure 8. Estimation of GPP and daytime Re components for the Yenicaga peatland based 
on Equation (6) (n = 1021 for NEE; and n = 896 for GPP and R E ). 
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Figure 9. Estimation of GPP and daytime Re components for the Yenicaga peatland based 
on Equation (1 1) (n = 1021 for NEE; and n = 896 for GPP and R E ). 




2 

A comparison of daytime Re versus ST max resulted in R a dj values of 39.3% and 99.3% based on 
Equations (6) and (11), respectively (Figures 10 and 11). 

Figure 10. A comparison of ST max versus daytime Re based on Equation (6) for the 
Yenicaga peatland (R 2 ad] = 39.3%; SE = 0.245; n = 896; P < 0.001). 
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Figure 11. A comparison of ST max versus daytime Re based on Equation (11) for the 
Yenicaga peatland {R 2 ad] = 99.3%; SE = 0.0039; n = 896; P < 0.001). 
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Soil water content was also determined to play a significant role in controlling daytime Re from the 
Yenicaga peatland. As a function of SWC, quadratic regression models elucidated 39.2% and 68.7% of 
variation in daytime Re according to Equations (6) and (11), respectively (Figures 12 and 13). A 
significant relationship between estimated GPP and PAR values was found, yielding R a( jj values 
of 80.3% and 55% from Equations (6) and (1 1), respectively (Figures 14 and 15). The rate of NEE in 

2 — i 

the Yenicaga peatland for the study period was on average estimated at 0.11 mg CO2 m s 
(26.5 kg C ha" 1 day" 1 ) and ranged from -1.26 mg C0 2 m" 2 s" 1 (-297 kg C ha" 1 day" 1 ) 
to 1.06 mg CO2 m" 2 s" 1 (249 kg C ha" 1 day" 1 ). This case points to the domination of total (positive 
daytime ± nighttime fluxes) Re over GPP (negative flux) throughout most of the study period, 
especially, towards the end of the period during which GPP decreased at a faster rate than i? E 
(Figures 8 and 9). During the peak growing season in July, GPP peaked and decreased as the season 
progressed. Over the study period, GPP was on average -0.74 ± 0.43 mg CO2 m s 
(-175 ± 102 kg C ha" 1 day" 1 ) and -0.35 ± 0.18 mg C0 2 m" 2 s" 1 (-2 ± 41 kg C ha" 1 day" 1 ) 
according to Equations (6) and (11), respectively. On average, nighttime R E was estimated 
at 0.23 ± 0.09 mg C0 2 m" 2 s" 1 (53 ± 22 kg C ha" 1 day" 1 ) based on the EC sensors, while daytime R E 
values were 0.64 ± 0.31 mg C0 2 m" 2 s" 1 (150 ± 74 kg C ha" 1 day" 1 ) and 0.24 ± 0.05 mg C0 2 m" 2 s" 1 
(57 ± 1 1 kg C ha" 1 day" 1 ) based on Equations (6) and (11), respectively. 
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Figure 12. A comparison of SWC versus daytime Re based on Equation (6) for the 
Yenicaga peatland (R 2 ad} = 39.2%; SE = 0.245; n = 896; P < 0.001). 
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Figure 13. A comparison of SWC versus daytime Re based on Equation (11) for the 
Yenicaga peatland (i? 2 ad j = 68.7%; SE = 0.025; n = 896; P < 0.001). 
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Figure 14. Quantification of relationship between PAR and GPP based on Equation (6) for 
the Yenicaga peatland (i? 2 adj = 80.3%; SE = 0.1929; n = 896; P < 0.001). 
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Figure 15. Quantification of relationship between PAR and GPP based on Equation (11) 
for the Yenicaga peatland (i? 2 ad j = 55.0%; SE = 0.1182; n = 896; P < 0.001). 
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Peatlands (bog hummock, bog hallow, and poor fen) in a cool temperate region of eastern Canada 
were reported to have NEE range from -0.18 to 0.13 mg CO2 m s ; nighttime R E range from 0.07 
to 0.36 mg C0 2 m" 2 s" 1 ; and GPP range from -0.20 to -0.33 mg C0 2 m" 2 s" 1 [15,16]. The reported 
ranges of NEE, Re and GPP values are in a close agreement with our findings for the Yenicaga 
peatland [NEE = 0.11 mg CO2 m~ 2 s _1 based on the EC data; nighttime Re = 0.23 mg CO2 m~ 2 s -1 
based on the EC data; daytime R E = 0.24 mg CO2 m s based on Equation (11); daytime 

2 —i 9 — 1 

Re = 0.64 mg CO2 m s based on Equation (6); GPP = -0.35 mg CO2 m s based on Equation (11); 
and GPP = -0.74 mg CO2 m s based on Equation (6)] , in particular, based on Equation (11). 
Similarly, EC measurements from 12 wetlands ranging from ombro trophic and minero trophic 
peatlands to wet tundra ecosystems across Europe and North America under temperate -to-arctic 
climate regimes showed that CO2 fluxes in July varied considerably between -3 and 1 g C m day 
for NEE; 1 and 4 g C m day for R E ; and -1 and -6 g C m day for GPP [17], which were close to 
our results. 

4. Conclusions 

Carbon metabolism components of the Yenicaga peatland as measured by EC sensors and 
quantified by M(N)LR models clearly revealed that diurnal and seasonal variations in exchange rates 
of CO2 between the atmosphere and the peatland ecosystem are dependent on the magnitude, rate and 
timing of GPP and total R E , which are in turn strongly controlled by the dynamics of soil moisture and 
temperature, and PAR. Peatlands experiencing drought conditions are reported to be able to act as a 
CO2 source. Given the total PPT: ET ratio of 0.3 during the study period, the Yenicaga peatland is 
considered to undergo a dry season and signals what future rate and amount of changes may be 
expected in the face of an increase in ET, air temperature and ecosystem disturbances as well as a 
decrease in PPT. The present study is the first one to quantify CO2 exchanges for a peatland in Turkey 
using real-time monitoring by the EC sensors. It is also the first time in Turkey that diurnal source and 
sink behaviors of a rarely occurring ecosystem like the Yenicaga peatland, which at the same time 
undergoes severe human-induced pressures such as conversions to rangeland and cropland, and peat 
mining, were quantitatively assessed using EC sensors. Further research is needed to explore an 
integration of remote and proximal sensors, and biogeochemical process-based models in improving 
our understanding and predicting impacts of human- induced disturbances on ecosystem metabolism. 

Acknowledgements 

We would like to thank TUBITAK and COST Action ES0903 for funding our project 
(CAYDAG109Y186). We are grateful to Mr. Aydin Giilmez and Staff of Yenicaga Municipality for 
their logistic support during our field studies. We acknowledge the financial and administrative 
contributions of Scientific Research Project offices of Abant Izzet Baysal and Akdeniz Universities for 
the project. 



Sensors 2011, 11 



537 



References 

1. Gorham, E. Northern peatlands: role in the carbon balance and probable responses to climatic 
warming. Ecol. Appl. 1991, 1, 182-195. 

2. Schindler, D.W. The mysterious missing sink. Nature 1999, 398, 105-106. 

3. Cayci, G.; Ataman, Y.; Unver, I.; Munsuz, N. Distribution and horticultural values of the peats in 
Anatolia. ActaHortic. 1988, 238, 189-196. 

4. Malone, S.; O'Connell, C. Ireland's Peatland Conservation Action Plan 2020: Halting the Loss 
of Peatland Biodiversity; Irish Peatland Conservation Council: Kildare, UK, 2009. 

5. Dengiz, O.; Ozaytekin, H.; Cayci, G.; Baran, A.; Characteristics, genesis and classification of a 
basin peat soil under negative human impact in Turkey. Environ. Geol. 2009, 56, 1057-1063. 

6. Sumer, N. Flora of Yeniqaga. M.Sc. Thesis, Graduate School of Abant izzet Baysal University, 
Department of Biology, Bolu, Turkey, 2002. 

7. Daughtry, C.S.T.; Gallo, K.P.; Goward, S.N.; Prince, S.D.; Kustas, W.P. Spectral estimates of 
absorbed radiation and phytomass production in corn and soybean canopies. Remote Sens. 
Environ. 1992, 39, 141-152. 

8. Running, S.W.; Hunt, E.R. Generalization of a forest ecosystem model. In Scaling Physiological 
Processes: Leaf to Globe; Ehleringer, J.R., Fields, C.B., Eds.; Academic Press: New York, NY, 
USA, 1993; pp. 141-157. 

9. Melillo, J.M.; McGuire, A.D.; Kicklighter, D.W.; Moore, B., Ill; Vorosmarty, C.J.; Schloss, A.L. 
Global climate change. Nature 1993, 363, 234-240. 

10. Thomas, M.V.; Malhi, Y.; Fenn, K.M.; Fisher, J.B.; Morecroft, M.D.; Lloyd, C.R.; Taylor, M.E.; 
McNeil, D.D. Carbon dioxide fluxes over an ancient broadleaved deciduous woodland in southern 
England. Biogeosciences 2010, 7, 3765-3814. 

11. Arneth, A.; Veenendaal, E.M.; Best, C; Timmermans, W.; Kolle, O.; Montagnani, L.; 
Shibistova, O. Water use strategies and ecosystem- atmosphere exchange of CO2 in two highly 
seasonal environments. Biogeosciences 2006, 3, 421-437. 

12. Billesbach, D.P.; Fischer, M.L.; Torn, M.S.; Berry, JA. A portable eddy covariance system for 
the measurement of ecosystem — atmosphere exchange of CO2, water vapor, and energy. J. Atmos. 
Oceanic. Technol. 2004, 21, 639-650. 

13. Coleman, K.; Jenkinson, D.S. ROTHC-26.3: A Model for the Turnover of Carbon in Soil. Model 
Description and Windows Users Guide; Rothamsted Research: Harpenden, UK, 2005; p. 45. 

14. Parton, B.; Ojima, D.; Del Grosso, S.; Keough, C. CENTURY Tutorial. Supplement to CENTURY 
User's Manual; Great Plain System. Research Unit Technical Report No. 4 [R]; NREC Pub. 
Natural Resource Ecology Laboratory, Colorado State University: Fort Collins, CO, USA, 2001. 

15. Moore, T.R.; Bubier, J.L.; Frolking, S.E.; Lafleur, P.M.; Roulet, N.T. Plant biomass and 
production and CO2 exchange in an ombrotrophic bog. J. Ecol. 2002, 90, 25-36. 

16. Lafleur, P.M.; Roulet, N.T.; Bubier, J.L.; Frolking, S.; Moore, T.R. Interannual variability in the 
peatland-atmosphere carbon dioxide exchange at an ombrotrophic bog. Global Biogeochem. Cy. 
2003, 17, 1036. 



Sensors 2011, 11 



538 



17. Lund, M.; Lafleur, P.M.; Roulet, N.T.; Lindroth, A.; Christensen, T.R.; Aurela, M.; 
Chojnicki, B.H.; Flanagan, L.B.; Humphreys, E.R.; Laurila, T.; Oechel, W.C.; Olejnik, J.; 
Rinne, J.; Schubert, P.; Nilsson, M.B. Variability in exchange of CO2 across 12 northern peatland 
and tundra sites. Glob. Change Biol. 2010, 16, 2436-2448. 

© 2011 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3 .0/) . 



